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Recently we have developed a theory of Keldysh formalism for mesoscopic systems. The resulting 
nonequihbrium Kubo formula for differential conductance makes it possible to propose the new 
formula of shot noise Sh, and thus to study shot noise in correlated systems at any temperature and 
any bias voltage. Employing this new approach, we analyze shot noise in the Kondo regime through 
a quantum dot for the symmetric case at zero temperature. Using the renormalized perturbation 
theory, we prove that in the leading order of bias voltage Sh equal to noise power at zero temperature 
conventionally used as shot noise. With Sh, we calculate the Fano factor for a backscattering 
current /;,: Fb = Shf^elt- It is shown that the Fano factor takes the universal form of Ft = 
1 + 4{R — 1)^/(1 + 5(-R — 1)^) determined by the Wilson ratio R for arbitrary strength of the 
Coulomb interaction. Using the Wilson ratio 7? = 2, our result coincides with the fractional value 
of Fb = 5/3 already derived in the Kondo regime. 



PACS numbers: 73.63.Kv 



INTRODUCTION 



Mesoscopic devices have made it possible to investi- 
gate many-body effect in transport phenomena. A bias 
voltage applied to a conductor drives the system into a 
nonequilibrium steady state. These systems thus offer 
promising opportunities to study the correlated trans- 
port out of equilibrium. 

Shot noise has been introduced to explore nonequilib- 
rium nature in these systems^'^:^. Originally shot noise 
was given by the zero temperature value of noise power S 
defined by the current-current correlation function. It is 
well known that, if uncorrelated particles with charge q 
rarely transmit through a conductor, shot noise is given 
by the Poisson value S ~ 2ql where / is the current. 
Therefore in this case, the quantity 5/2/ has a possibil- 
ity to determine the unit of charge carriers q, extracted 
from not only the current^'^'^. In literature, this result 
is specific to a low transparent conductor without cor- 
relation effects. However, the possibility of shot noise 
has stimulated recent studies in other systems, for in- 
stance, correlated systems. There, the ratio S/2I has 
been estimated as an indication of effective " charge" . As 
a successful example to give a true quasiparticle charge, 
the fractional charge e* = e/3 in the fractional quantum 
Hall regime was determined through S/21'^'^ . 

Shot noise defined by the noise power has been useful 
to reveal attractive features at low temperatures. Exten- 
sion of the concept of the shot noise at finite tempera- 
tures may open new possibilities out, however, how do 
we define shot noise at finite temperatures? 

To illustrate the problem we take a quantum dot 
system as a stringent example^. In the noninteract- 
ing case, noise power splits into thermal noise St{T'^) 
and shot noise Ss{T°) as 5° = StiT°) + S,{T°) at any 
temperatureii^i^, where T^ is the bare transmission prob- 
ability. The fact enables us to define the shot noise 



as S'^ — St{T^) at any temperature and any bias volt- 
age. With correlation effects, noise power is obtained as 
S = StiT) + Ss{T) + AS. If the effect of Coulomb in- 
teraction may be included by simply changing T° to the 
full transmission probability T, noise power would be ex- 
pressed as St{T) + SsiT). Shot noise thus could be de- 
fined in the same way as the uncorrelated system. How- 
ever, an additional term of AS* contributes to S. From 
the expression of AS*, it is not possible to split it into 
the thermal- noise part and the shot-noise part. Thus in 
a quantum dot system it has not been possible to extract 
shot noise from the noise power. Until this point, we have 
stated the problem for a quantum dot system. However, 
the problem to define shot noise in correlated systems is 
a general one, not limited to the quantum dot system. 

Recently we have succeeded in extending a theory of 
Keldysh formalism for mesoscopic systems, and thus pro- 
vided a basis for the definition of shot noise at any tem- 
perature. 

In this study, we have focused on a group of mesoscopic 
systems, typically designed as a correlated conductor at- 
tached to extended reservoirs through tunnel barriers. 
We have specifically examined the density matrix of the 
Keldysh formalismiiSiS for these systemsiSiliii^. An at- 
tempt to obtain an explicit expression of the density ma- 
trix had been already made in Refs llSlllil . However, this 
attempt was not fully successful, since the density ma- 
trix still contained infinite series of operators which arose 
from the formal expansion of the s-matrix. Summing up 
them, we have thus proven that the density matrix has a 
form proposed by Maclennan and Zubare v^^i^^'^^ . This 
type of density matrix was also derived by using a C* 
algebrapiSiiS. 

The resulting density matrix allows us to derive a gen- 
eral expression of differential conductance: G = PS/A — 
(3Sh/4:, where S is the current-current correlation func- 
tion, and Sh is the non-trivial current-charge correlation 



function and (3 — l/k-oT^'' . We will call it the nonequi- 
librium Kubo formula. Furthermore, the density matrix 
has made it possible to prove that a steady state is re- 
alized, and describe the issue of irreversible processes in 
Keldysh formalism. Therefore our theory gives a unified 
description of desirable features in mesoscopic systems 
out of equilibrium. 

Moreover, the nonequilibrium Kubo formula allows us 
to introduce a natural definition of the shot noise in 
correlated systems. In fact the current-current corre- 
lation function S is noise power, and of course G is 
differential conductance. With these observable quan- 
tities, the nonequilibrium Kubo formula is written into 
Sh ~ S — ikBTG, and thus we propose the current-charge 
correlation function Sh as the new definition of the shot 
noise. This concept has been confirmed in a noninteract- 
ing quantum dot by explicitly calculating S, G and the 
new formula Sh with the Keldysh Green function. The 
resulting relation among physical quantities is called as 
the nonequilibrium identity. Thus we can compare the- 
oretically calculated Sh at a certain temperature with 
S — AksTG determined from S and G measured in ex- 
periments. Therefore, the nonequilibrium Kubo formula 
makes it possible to give a new perspective on studies of 
shot noise at any temperature in correlated systems»ii. 

Our aim in this paper is to apply this new approach 
to the quantum dot system in the Kondo regime by an- 
alyzing shot noise at low bias voltages. 

First, we briefly sketch recent studies of shot noise in 
the Kondo regime. In these studies shot noise is defined 
by the noise power S at zero temperature. Theoretically 
shot noise was investigated in the s-d limit^''. When a 
bias voltage eV is larger than the Kondo temperature 
Tk, with decreasing the bias voltage shot noise increases 
logarithmically as is well known in the Kondo systems. 
Close to the unitarity limit of eV < Tk, shot noise is 
suppressed as 






h 3' 'ITk 



(1) 



where T^ji are the resonance widths between dot and 
left and right leads. The Fermi-liquid fixed point Hamil- 
tonian H—— is used in calculations. 

In recent studie a^'^i^^ , the Fermi-liquid features in 
eq.([T]) have been reexamined at zero temperature. The 
fixed point Hamiltonian is characterized by a one-particle 
scattering and a two-particle scattering. They give 
corrections to the 7r/2 phase shift in the unitarity 
limit: a frequency- liner phase shift and a quasiparticle- 
distribution- linear phase shift. 

As a result, a current h = —edt{NLH {t) — Nrh {t)) de- 
scribes a reduction from the perfect transmission for the 
7r/2 phase shift in the unitarity limit, which is the defi- 
nition of the backscattering current in Refs ]23ll24l . Here 
Nl,r express particle-number operators in left and right 
leads and A^it) is the Heisenberg representation for an 
operator A. 



At zero temperature, the observable current is defined 
by / = 1e}V jh — If,, where /„ — 2e^V/h is the current 
in the unitarity limit. Since 2e^V/h becomes a constant 
value, noise power given by the current-current correla- 
tion function is rewritten as 



f 



dt{{5hH{t).5hHm)- 



(2) 



Therefore the shot noise is determined by only the 
backscattering current when eT^ < Tk- 

On the basis of this fact, using S/2Ib, it has been 
attempted to estimate an effective "charge" for the 
backscattering current^^. In the present work we will 
discuss the Fano factor Ff, = S/2eIb in units of charge e 

In Refl23l. the symmetric case F^ = Tr is especially 
focused on. This condition makes the leading order of S 
in the bias voltage 0{V^) as in eq.((T]) . Then, the leading 
order of Ii, also becomes 0{V^). In the symmetric case, 
Fi, — S/2eIfj results in a universal fractional value up to 
0{V^) 



Fh = 



S 
2^b 



(3) 



It has been pointed out that this universal enhancement 
from the unity originates from the two-particle backscat- 
tering. 

This universal feature has stimulated further studies: 
in the context of the full counting statistics^^ this result 
has been reproduced, and a shot-noise measurement has 
been reported on F^ = 5/3 close to the unitarity limit at 
low bias voltages^. We address this topic from our point 
of view. 

In this paper, we begin with a brief review of the gen- 
eralization of the Kubo formula into a nonequilibrium 
situation and shot noise in SecO In Sec|3] applying our 
new concept to a quantum dot system described by the 
Anderson model, we investigate the shot noise and the 
Fano factor in the Kondo regime. We concentrate on the 
symmetric case at zero temperature. At zero temper- 
ature the nonequilibrium Kubo formula yields S — Sh- 
Thus, using the renormalized perturbation theory, we cal- 
culate noise power S and shot noise Sh independently up 
to the leading order for the symmetric case: 0{V^) and 
show that they are identical. The backscattering current 
is also calculated up to 0{V^). Finally we will discuss 
the universal feature of the Fano factor for any Coulomb 
interaction up to 0{V'^). 



SHOT NOISE AND NONEQUILIBRIUM 
KUBO FORMULA 



We start with a brief summary of shot noise and the 
nonequilibrium Kubo formula based on Keldysh formal- 
ism for mesoscopic systems in nonequilibrium steady 

states-!^. 



2.1. density matrix 

Our starting point is the Hamiltonian for a class of 
mesoscopic systems: a correlated conductor He attached 
to left and right infinite reservoirs Hlr through bound- 
ary couplings HcL.R- For instance a quantum dot system 
is a typical example of them. The Hamiltonian is de- 
scribed as 



H = Hc + Hl + Hr + H, 



cLM- 



(4) 



He is given by noninteracting and interacting terms as 
He = Hco + Hci- Each reservoir must be sufficiently 
large to behave as a good thermal bath. 

In a conventional explanation of Keldysh formalism, it 
seems that concepts such as the Keldysh contour, and the 
Keldysh Green function are stressed^. In contrast, our 
main argument is to examine the early stage of Keldysh 
formalism where the perturbation term is adiabatically 
turned on. 

To begin with, we divide H into an unperturbed term 
Hq = Hco + Hl + Hji and a perturbation term Hi = 
HcL,n. + i?ci— '^'^- To adiabatically introduce the per- 
turbation term, we use a time-dependent Hamiltonian as 
H,t = Ho + e-'\'\Hi. 

The initial state at io = ~oo is determined by the three 
separated systems in Hq'. left and right reservoirs and the 
conductor which have different chemical potentials ^il.Ri 
^c- We define a bias voltage as eVl2 = (^l — A^fl)/2 
and an averaged chemical potential /x = (/i^ -I- /iij;)/2, 
where /i = /^c is assumed. The origin of energy is set to 
fi = fic = 0- As a result the initial density matrix is given 
by 



Po = e 



-0{H„-eV/2{NL-NR)-no) 



(5) 



[Ho,Nl — Nil] = is assumed throughout. The initial 
thermodynamic potential is chosen to satisfy Tr{/9o} — 1, 
no = -l//31nTre-'3(^o--^/2(^--^H)). 

The time evolution of the density matrix p^ (t) is deter- 
mined by the Neumann equation. With Pe{t), the expec- 
tation value of any operator O is defined by Tr{p^{t)0}. 
It is transformed as 

Tr{p,(i)O} = Tr{p,(0)O^(i)}, (6) 

where the Heisenberg representation is given by Oh (t) — 
Se{Q,t)O{t)S,{t,0) and 0{t) = e'^o*C'e~'^°*. 

Here in order to avoid confusion, we comment on the 
time dependence of the expectation value. It is deter- 
mined by only the Heisenberg representation of an opera- 
tor O as shown in the right-hand side of eq. (O . Concern- 
ing the density matrix, it is sufficient to consider pe{0) 
at t = 0. We thus consider a specific time of i = 0, but 
until this point we do not use any assumption either the 
steady state is realized at i = or not. After taking a 
limit of e ^ 0, the system reaches a steady state. Later 
we will return back to this point in Sec l2.4l 



We thus focus on the density matrix Pe(0) = Pe which 
is given by a formal solution of the Neumann equation 

Pe = Se{0,-Oo)poSe{-CO,0). (7) 

Using the unitarity of Se{0,—co)Se{—oo,0) — 1, p^ be- 
comes 

p, = exp{-/3(5,(0,-oo)ilo^e(-oo,0) 

-VS,{0, -oo)e/2(7Vi - NR)S,i-^, 0) - r!o)}. (8) 

To proceed calculations, we have to compute a type of 
quantity 



A, 



S,{0,-oo)ASei-oo,0). 



(9) 



We have two cases in our mind: A = Hq and A = 
e/2{NL — Nr), where we suppose that A satisfies 
[A, Hq] = 0. The critical step in calculations is to deriveii 



A, = A+ / dte-'^'^JAHit). 

..' — oo 

The "current" for the operator A is defined as 

d 



JAH{t) 



dt 



Auit). 



(10) 



(11) 



Technically eg. pU)) is sufficient to calculate the den- 
sity matrix. However, we shortly discuss eg. pTI)) . leading 
to deeper understanding of the nonequilibrium nature in 
Keldysh formalism. Integration by parts in eq. (|10|) yields 



A,^e dte-^l*Uff(i). 

J — oo 



(12) 



Thus A^ expresses a long-time average of A in the limit 
e — > 0. In fact, A^ defined by eq. (|12p represents noth- 
ing but the invariant part of an operator A^ introduced 
by Zubare\»i^. We have thus proven that the adiabatic 
introduction of Hi in Keldysh formalism corresponds to 
taking the invariant part by Zubarev. We revisit this con- 
cept to discuss steady states and irreversible processes in 
Secim 

We use eq. pH)) in two cases: A = Hq and A = 
e/2{NL — Nji) where the "currents" are defined by energy 
change Jeuit) and charge current JcH{t) using eq. pTj) . 
respectively. Substituting them into eq.® yields 



exp 



.0 
l3{Ho+ dte-^l*! JeH(t) 

J — OO 

eV/2{NL -Nr)-V dte-^l*l Jeif (t) - n^ 



(13) 



The density matrix of Keldysh formalism thus becomes a 
type of the nonequilibrium statistical operators initiated 
by MacLennan and Zubare v^^i^^ . This form was also 
obtained by the C* algebr a^^i^^ . 



2.2. nonequilibrium Kubo formula 

The resultant density matrix makes it possible to gen- 
eralize the Kubo formula for conductance into a nonequi- 
librium steady state. For convenience, in this work we 
introduce a bra-ket notation: (O) = limTrlpjO}. The 

current is thus denoted as 



e^O 



/^(J«(t)) = (J), 

where in the last equality of eq. (|14p . 
ture later proved in Sec l2.4l is used, 
differential conductance is given as 



(14) 

a steady-state fea- 
From the current, 



(15) 



Eq. p4|l shows that the current / depend on a bias voltage 
V only through the density matrix p^. To obtain G, 
dpc/dV should be calculated. Thus G is given by 



G = lim/3Tr<j j( e/2{NL - Nr) 



/ dte-^l*ljH(t) 

•J — OO 



drtp 



(16) 



Using the initial thermodynamic potential fJo and 
S^Sl — 1, dVlp/dV is calculated as 



lL«sf? = -/^(»' 



Nr. 



dte-^l*l(JH(t)). 



Thus the inside of (...) in eq. (fT6|) is expressed by a fluc- 
tuation of an operator: 5Af^ = A^ ~ i^e)- In fact, A^ 
just expresses the invariant part oi A = e/2(7Vi — Nr) in 
eq. (fTO|) . For simplicity we abbreviate it A^ for a while. 
Furthermore, using (BSG) = {6B6C) for any operators 
B and C, ea. p^ is rewritten into G = (3{SJSA^) 

Here we symmetrize the expression of G = (3{5J5Af). 
The current operator J and A,, do not commute, but it 
is possible to transform to the inverse order of {5A^5J). 
Using S'e[i?o, e/2{NL - Nr)]SI = 0, A^ is proven to com- 
mute with pe. With a cyclic permutation inside the trace, 
G = (3{5J5A^) can be written into G = (3{5A^5J). We 
have thus obtain two different expression of G. 

Moreover each form of G is characterized by the time- 
integral over negative t. It is rewritten into one over 
positive t using the steady-state feature which will be 
discussed in eq. p4| . Consequently, four different expres- 
sions of G are derived. Symmetrizing them yields 

G^^{{5J,e{5NL~5NR)}) 



where {B^G} 
tion relation. 



di({<5J^(t),<5J^(0)}), (17) 



BC + GB represents the anticommuta- 



We have found that differential conductance G is de- 
termined by a current-current correlation function, and 
an unusual current-charge correlation function. In the 
linear response regime for G: ^ = 0, we can show that 



({(5J,e((5iVL-<5A^K)})v=o = 0, 



(18) 



and thus eq. lfTTl) reduces to the standard Kubo for- 
mula. Away from the linear response regime, the current- 
current correlation function gives a naive extension from 
the relation between fluctuation and dissipation. On the 
other hand, the current-charge correlation function gives 
an intriguing modiflcation far from the linear response 
regime. In conclusion we have succeeded in obtaining a 
nonequilibrium Kubo formula in mesoscopic systems. 



2.3. shot noise 

We reexamine the nonequilibrium Kubo formula in 
view of physical quantities. As measurable quantities, 
there are differential conductance G and noise power 
given by the current-current correlation function S 



S^ 



dt{{5JH{t),6JH{0)}). 



(19) 



At a steady state noise power at zero frequency domi- 
nates. 

On the other hand, we define the current-charge cor- 
relation function as Sh 



Sh = -{{SJ,eiSNL-dNR)}). 



(20) 



The nonequilibrium Kubo formula in eq. (|17p leads to 
the relation among Sh and two observable quantities G 
and S 



Sh = S- AkBTG. 



(21) 



As stated in the introduction, in uncorrelated systems 
shot noise is given by S° — St{T°)- Considering the fact 
that the thermal noise 5't(T") = 4/cbTG°, Sh in eq. ipT]) 
suggests the natural extension of the shot noise. There- 
fore, we propose to define the shot noise by the current- 
charge correlation function Sh at arbitrary temperature 
and bias voltage in correlated systems. 

As a result, eq. (PT|) expresses the relation among phys- 
ical quantities out of equilibrium. Thus we will call it a 
nonequilibrium identity. In the linear response regime, as 
discussed previously the current-charge correlation func- 
tion vanishes, and thus Sh ~ 0. Eg. ipT]) then goes back 
to the Nyquist- Johnson relation S{0) = 4fcBTG(0). On 
the other hand, at zero temperature eq. (j2ip shows that 
Sh indeed equals to S* at T = which is nothing but the 
original definition of the shot noise. 



2.4. steady state 

Having discussed consequence of the density matrix 
on the transport property: the nonequihbrium Kubo for- 
mula and shot noise, here we turn to others: irreversible 
processes and steady states in Keldysh formalism. 

We have shown that eq. pI7|) . more explicitly eq. fT^ 
corresponds to the invariant part of an operator A dis- 
cussed by Zubarev. Originally, this concept was intro- 
duced in a local equilibrium system. Zubarev could show 
that the entropy production was positive in the linear 
response regime, and thus concluded that the method 
could describe dissipations. Following the same analysis, 
it is possible to prove that Keldysh formalism have the 
same feature as Zubarev theory. 

Secondly we discuss steady states. We have shown in 
RefllTI that the commutation relation between the invari- 
ant part and the Hamiltonian satisfies 



\im[Ae,Het] 

£—►0 



lim —ieiAf 



A). 



(22) 



Assuming that A^ exists in the limit e ^ 0, the r.h.s. of 
ea. (l22p vanishes. In e — > the invariant part commutes 
with the Hamiltonian, and thus conserves. 

This result is crucial for the proof that a steady state 
is realized in the limit e -^- 0. Using a similar technique 
to derive eq. pUj) , the expectation value of an operator O 
is rewritten into 

lim Ti{p,0 Hit)} ^ limTrjOeC} 

+ lim / dt Tr{i[p„ H^t']JoH{t - t )}. (23) 

The definition of p^ in eq. ([7]) shows that the density ma- 
trix itself in fact becomes an invariant part. As a con- 
sequence the commutation relation \p^^}i^^i\ in eq, 
vanishes in the limit e ^ 0, leading to 



limTr{peC//(i)} = limTrjp.O}. 



(24) 



We have proven that a steady state is realized in the limit 
e — > 0. According to the same analysis, (C'_H'(t)0_f/(i')) — 
{Onit - i')Cff (0)) can be also derived. 



3. SHOT NOISE IN THE KONDO EFFECT AT A 
QUANTUM DOT 

3.1. Quantum dot system and physical quantities 

Now let us turn to the discussion on shot noise through 
a quantum dot in the Kondo regime. We consider the 
Anderson model with a single-level, 

koL(7 <y 

kotcr 



dl and c{.Q,^(a = L,R) create an electron with spin a at 
the dot and the left-right lead respectively. The last term 
describes the tunneling between the dot and the leads, 
which determines the resonance width T(uj) = {Tl(uj) + 
rfl(a;))/2 withri,ii(tj) = 2TrJ2k |HL,fl|''5(w-efcL,fl). In 
the limit of large band width, the resonance width may 
be assumed to be a constant F = (F^ -I- Tfj)/2. 

As discussed in the introduction, we focus on the 
Kondo effect near the unitarity limit at zero tempera- 
ture, and investigate shot noise up to 0{V^). To discuss 
the unitarity limit, it is natural to concentrate on the 
symmetric case, e^ = —U/2 and F^, = F/j = F. 

Traditionally noise power S at zero temperature has 
been defined as shot noise. In contrast we define the 
shot noise in general by eq. (|20p . From the nonequihb- 
rium identity eq. (|2ip based on the nonequihbrium Kubo 
formula, 



Sh = S, 



(26) 



holds at zero temperature. Our first aim in this paper is 
to show that S = Sh holds up to C'(F'^) by explicit calcu- 
lations of the two quantities. Second aim is to calculate 
the Fano factor^i^iiSS^ 



Fh = 



Sh 
2e/h' 



(27) 



where the backscattering current lb is defined as It = 
lu ~ I with the current / and the current in the unitarity 
limit /,, . 



3.2. Renormalized perturbation theory 

For explicit calculations we employ the renormal- 
ized perturbation theory (RPT). It was introduced by 
HewsoriSIi^s j^ equilibrium, and then extended by Oguri 
to the quantum dot system under a finite bias voltage—. 
In this section the RPT applied for the Anderson Hamil- 
tonian in eq.(|25p is introduced first in equilibrium, and 
then out of equilibrium. 

The perturbation theory in U for the Anderson model 
is characterized by a set of parameters: the energy level 
Cd, the resonance width F, and a regularized Coulomb 
interaction U/ttT. In the RPT these parameters 
are substituted by renormalized ones, e^, F, U/ttT — > 
id, f , tJ /nf. Using ward identities, three parameters are 
proved to be related with the regularized spin suscepti- 
bihty Xs = [2TrT/{gpB)'^]Xs, charge susceptibihty Xc = 
(7rF/2)xc and specific heat coefficient 7 = (3F/27rfcg)7, 

Xs = p.(0)(l + [7p(0))7rF, 

Xc = p(0)(l-C7p(0))^F, 

7 = P(0)^r. (28) 

The density of states at the Fermi level is given by 



ttF (e<i/F)2 + 1 



(29) 



Eq. ([28| leads to the Fermi liquid relation: Xs + Xc = 27, 
thus two quantities Xs and Xc become independent vari- 
ables. The Wilson ratio defined as i? = Xs/l is expressed 
as -R = 2xs/{Xs + Xc)- As another relation, we consider 
the Friedel sum rule 



1 1 
tan 

2 n 



■© 



(30) 



Thus Xs, Xc and n determine the three renormalized pa- 
rameters 



f 

U/ttT 



-sin(2n7r)r/(xs + Xc), 
2sm\nTT)T/{x,s+Xc), 
{R- l)/sin2(n7r). 



(31) 



From the bare Hamiltonian the renormalized parameters 
can be calculated using the exact Bethe ansatz results of 
Xs, Xc and f»^. Alternatively these parameters can be 
also estimated by the numerical renormalization group^. 

Here we comment on the behaviors of these parameters 
in the s-d limit and in the symmetric case. 

In the Kondo regime, according to Nozieres^, we de- 
fine the Kondo temperature as Tk = (g^J.B)'^ /t^Xs- The 
Bethe ansatz result of Xs allows us to determine the ex- 
pression of the Kondo temperature^! 



rK^iyZ7fcxp(-!^*i±^ 



(32) 



In the Kondo limit, n — > 1/2 and Xc ~> 0, therefore 
R = 2, Ed = 0, and 7rf = U. Moreover f = 2r/xs = 
[g^jL-a)^ /t^Xs is derived. Using the definition of the Kondo 
temperature, F = Tk is obtained. 

The symmetric case gives rise to n = 1/2, thus e^ = 0, 
and 



^ — '■ 



(33) 



Consequently physical quantities are expanded in a 
power series of i? — 1 for the symmetric case. 

As a check on RPT, the second-order correction in U 
to the self-energy has been calculated up to 0{uP') and 
0{T^) for the symmetric case^las 



'^^ ' 2r VttF 



{uj^ + ir^T^) + 



(34) 



This result is in agreement with the one in ReflSlI where 
all orders in the bare U are calculated at low frequencies 
and temperatures. Within 0{uj'^) and 0{T'^), the second- 
order contribution in the renormalized U thus gives the 
exact expression. RPT enables us to discuss the exact 
Fermi-liquid features in the low energy region for all pa- 
rameter regimes of e^, F, and U/nT from the weak cor- 
relation regime to the strong correlation regime. 

Having discussed RPT in equilibrium, let us turn to 
extension into under a finite bias voltage^S. Using the 
renormalized parameters, this procedure is done as fol- 
lows. 



First we introduce Green functions based on Keldysh 
formalism 

Gd(T-0) = -i(T,d,(T)4(0)), 
GpUr - 0) = -i{TcCpkAT)dlm, 

Grfpfc(T-o)^-z(r,d.(T)4,(o)), 

Gpkp'k'ir - 0) = -i{TcCpka{T)c^p, 1^,^(0)), (35) 

where t = t is the time variable. The lower and upper 
branch along the Keldysh contour are denoted by — and 
+, and thus Keldysh components are defined as A(i" — 
0'3) = A"^(i-0) ioT a,(3 = -,+. 

The equation of motion enables us to relate these full 
Green functions with Gd 

Gd{u:) ^ 9°d{u:) + 9°dHnuo)Gd{uo) 
Gpkd{uj) = Vpkgpkpk{^)'^zGd{cj) 
Gdpki^) = Gdi^)<7zV*,,g°i,pki^) 
Gpkp'k'i^) = 9pkp'k'i^)^PP'^kk' 

+ Vpkgpkpki^)'^zGdiuj)a^V*,i,,g°,k,p,k'i^)- (36) 

G represents the matrix (G)"'^ = G"^ and (Tz is the third 
Pauli matrix in the a/3 space. Here g^ refers to U = 
and VL,Rk = 0. 

As the zeroth-order Green functions in U, we consider 
the initial Green functions g = GIp^Q, renormalized by 
VL,Rk- The initial Green functions are obtained as the 
solution of the Dyson eq. for Gdlfj^o = 9d in eq. 



9d'{^) = (1 - UsWdi^) + fcS9di^) 

ffrf+M = -/cff(3SH-3dM) 

-g+-{Lo)^il-UWdi^)-rd{u^)) 

-g++{uj) = -(1 - .UsMi^) - Ufd{^). (37) 



;(^) = 1/(0; -e-,-fzf),g^(^)=g- (.;);, 



/eff = 



F, 



(38) 



Here the Fermi distribution function Jl.r = 1/(1 + 
g/3(w=Fey/2)^ are not renormalized. Substituting Gd\u=o — 
gd into eq. (I36p , we get the explicit forms of other initial 
Green functions gpkd, gdpk and cjpkp'k'- In the pertur- 
bation theory in U, we use these initial Green functions 

9- 

Using the initial Green function gd, the Dyson equation 
of Gd in eq. (j36p is rewritten into 



Gd{Lj) = gd{Lu) + gd{uj)T.diuj)Gd{Lu). 



(39) 



With a unitary transformation, eg. ([55)1 is reduced to the 
three-components form. One of them becomes the Dyson 
equation for the retarded component. 



G5(c.)=55(c^)+5SMS^MGS(c^). 



(40) 



The current is expressed as the renormahzed form, 

/oc 
du;T{cj){h^f^), (41) 

-°° 

where transmission probability T(w) is defined as 



Tico) 



TrV 



Li R 



Tr 



■p{uj), p{uj) 



1 



ImG^(c.) (42) 



Under a finite bias voltage, the second-order calcula- 
tion in tj for 17 has been done for the symmetric case^S, 



E5(c.) = z 



1 

2r 



£_ 

nf 



[lo'^ + {eVf + tt'^T'^) + ■ ■ ■ .(43) 



Clearly this corresponds to eg. ([34]) in the equilibrium 
limit F = as expected, and moreover reproduces the 
complete result if combined with the Ward identity^. 

We substitute S''((^) of eg. ([43] ) to the formal solution of 
G^ in eq. (|40| , and calculate the transmission probability 
with eq. 



T(u;) 



1 



f2 



1 

2r2 



u 

7rf 



LO 



-Mvf 



(44) 



The resulting transmission probability allows us to cal- 
culate the current defined in eq. (|4ip . For later use, we 
show the current at zero temperature 




where U /ttT = i? — 1 is used in eq. ([45|) . In the Kondo 
limit, R — 2 and T = Tk- The s-d limit result given in 
eg . (j46l) completely agrees with the ones in Refs l22ll23ll24l . 
I25I where the fixed point Hamiltonian are used. 

Therefore, for the quantum dot system the second- 
order calculations in U to the self-energy can provide the 
exact result up to 0{uj^), ©(T^) and 0{V^). This result 
leads to the correct expression of current up to 0{V^), 
where we need not count higher-order corrections in U. 



3.3. Calculation of S and Sh 

In this section, we consider shot noise for the symmet- 
ric case at zero temperature up to 0(V^), by employing 
the two definitions of the shot noise, noise power 5* at 
T — 0, and the new formula for shot noise Sh- Within 
0{V^), in the same way as the current, the second-order 
calculations in U are sufficient to give correct results for 
S and Sh- 

Both the noise power S in eg. p^ and the shot noise 
Sh in eq. (PD|) are two-particle Green functions. Gen- 
erally they are expressed by the bubble-type diagrams 



(bubble-diagrams) and the vertex-correction-type dia- 
grams (vertex-diagrams). We use subscripts 1 and 2 
to denote contributions from the bubble-diagrams and 
vertex-diagrams. Thus S = Si+ S2 and Sh = Shi + Sh2- 



(a) 



(b) 









(c) 



(d) 







FIG. 1: For the symmetric case, the second-order corrections 
to the vertex-diagrams. 



3.3.1- Noise power S 

This section is devoted to the analysis of noise power 
S. As mentioned previously, RPT in JJ based on Keldysh 
formalism is employed. We begin with giving a definition 
of S in the Keldysh form. Considering the steady-state 
feature in eq. ([M)) . the original S in eq. p^ expressed by 
the anticommutation relation is rewritten into 



dt{ T, JH{t+)JH{0-))l 



{2Re{K+--K+-)} 



(47) 



{■..)l means the linked parts. Tc is the time-ordering op- 
erator defined on the Keldysh contour. The time variable 
t represent times on the upper and lower branches. It 
is convenient to classify the current-current correlation 
functions into K^'^ 

Kb=JdtK+^{t-0), (48) 

where if+^(i - 0) ee Ka,b{t+ - 0"). Here 

pp' kk'aa' 

X ( Tc C^HpkAr)dHAT) C^Hp'k'aMdHa'mL, 

pp' kk' aa' 

X ( Tc c^HpkAT)dH.{r) dL,(0)cffp.fev'(0))L 

Ka,b{T — 0) are also the two-particle Green functions. 
Thus Ka.b is generally given by contributions of the 
bubble-diagrams Kia^b and the vertex-diagrams K2a,b as 
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Ka^b = Kia.b + K2a,b- Appling this relation for S in 
ea. (|T7)) . we define S'1.2 to satisfy 5 = S'l + ^2 as follows, 



5i = ^{2Re{K+--K+-)} 



S2 



{2Re{K+--K+-)} 



(49) 



Thus, explicit calculations of Ki,2a,b determine noise 
power S. 

Let us start with the discussion of Si. Kia^ are eval- 
uated as 



i^ia(T-O) = (^V2) 



pp' 



/ ^ VpkVpik' 

kk' aa' 



X Gdp'k'{T~0)Gdpk{0-T)S, 



Kuir-O) = {f/2) 






aa' 1 
kk' aa' 



X GdiT-0)Gp.k'pk{0~T).S„^, (50) 

In the following, we use the Fourier representation. 

It is shown that Gdpk and Gp>k'pk are related with Gd 
in eq. (l36|) . Technically, these relations are sufficient to 
proceed further calculations. However, here it is better 
to comment on a difference between calculations of noise 
power and shot noise. 

Here we concern with the /c-sum in ea. (|50p . In both 
cases the fc-sum is reduced to a form 



E 



VpkGdpk , / VpkVp'k'Gp' 



k'pk 



kk' 



k 



\Vpk\ dpkpki^)- 



dpkpki'^) ^'^^ ^^6 bare Green functions for the free elec- 
trons in leads when [/ = and Vpk ~ 0. Thus the ana- 
lytic properties of 3°^, ^.(a;) are characterized by infinites- 
imal quantities ±ie. By performing the fc-sum of g^kpk^ 
the ^-function singularities are averaged. Therefore in 
the full bubble-diagram calculations, the fc-sum provides 
well-defined quantities. 

However, in the shot-noise calculations, another type 
of the fc-sum leads to a ^-function singularity. We will 
discuss this point in the next section. 

Let us return back to the main discussion. The rela- 
tions in eg. ipS)) allows us to rewrite 2 ^{^{K^' — K^^) 
into an expression with T, 

+ATiu;)ifLil-fR) + Ml-fL)). (51) 
Using 

hil - Ir) + fnil - h) = 

fL{l-fL) + fR{l-lR) + {fL-fR)\ (52) 

Si can be represented as 



^1 



4e^ 



dujf{u){fL{l-lL) + fR{l-fR)) 



4e^ 



+ ^ / dujT{u;){\ - T{u;))ifL - IrT- (53) 



We have proven that noise power Si for the bubble- 
diagrams decouples into the thermal-noise part and 
shot-noise part in the same manner as noninteracting 
systemsi"^. However, in this case the transmission prob- 
ability is fully renormalized by the Coulomb interaction. 
At zero temperature, in Si the thermal-noise part van- 
ishes. Consequently we discuss the remaining shot-noise 
part in the second- line of eg. ([55]) . With T given by 
ea. (|44)) . we evaluate the asymptotic form of Si 

where U /-kT ^ R-l. 

Here, let us turn to the discussion of S2 for the vertex- 
diagrams. Diagrams up to the second order in U are 
shown in Fig. [TJ For the symmetric case, we can show 
that the contribution from the diagram in (c) cancels out 
the one from (d), and the one from (b) itself vanishes. We 
explicitly calculate the remaining diagram in (a). 

We calculate K2a up to the second order in tJ 



K2a,b{T - 0) = 2 / dTidT2Fa,b. 



(55) 



where 



Fa = Mr - n) E ^z"' E Vp'k'gap'k' in - 0) 
p' fc' 

X gd{0-T2)^aPP^Vpkgdpk{T2-T) 
P fc 

X Ugd{Ti - T2)Ugd{T2 - Ti), (56) 

Fb = gd{T - Ti)gd{Ti - 0) X E 

pp' 

E K'k'9p'k'd{0 - T2)aP'P E ypkgpkd{T2 - t) 

k' k 

X Ugd{Ti - T2)Ugd{T-2 - Ti). (57) 

With eq. (|49| 5*2 is obtained in the Fourier representation 

4e2 /■ 
5*2 = — c / dujiduj2 Fab, 

where c — i7^f^/(27r)^. Fab is given by 

Fab = Y,gd{i^ir-<yT9d{^l)^''{lL - /fl)(^l +C^2) 
a 

X 5d(^2)"".gd(t^2)"+(./L - ./fi)(^2) X T{UJ2 - UJi), 

where a = -a and (/l - /fl)(w) = /l(c^) - /i?(cj). The 
vertex function r(a;) is defined as 



T{u:)^ j dw'g-+{u:')g+-(u: + u:') 



(58) 



Performing the frequency-integral gives the explicit form 
of r(a;). In eg. ([55]) we show the leading-order of r(a-'). 
Finally S2 becomes 



^-^iHi'-'Kt' 



(59) 



Therefore S — Si + S2 iS also obtained up to the 
leading-order in V, 



->U(¥)-I<-'K¥ 



liR-^)"('-f) >-m 



3.3.2. Shot noise Sh 

This section is devoted to calculations of the new for- 
mula of shot noise 



Sh^-{{5J,e{5NL~5Nn)}). 



(61) 



Originally, Sh contains 4x4= 16 two-particle Green 
functions. The current conservation ( J^ -I- Jp) = re- 
duces the number to 8 two-particle Green functions. As a 
preliminary step, we define the reduced Sh in the Keldysh 
form 

Sh = ^2 Re{{KLL + Krr) - {Klr + Krl)} , (62) 
where Kppi for p,p' — L, R are given by 

PP' ^ / ^ ^pk 

kk' aa' 

^{Tc 4p,,(0+)dff,(0+) C^Hp'k'.'iO')cHp'k'a'iO-))L. 

The equal-time correlation functions arc defined by using 
the Keldysh branches. 

Kppi as the two-particle Green functions are generally 
given by Kippi for the bubble-diagrams and K2pp' for 
the vertex-diagrams. Correspondingly, 5^1,2 are defined 
as Sh = Shi + Sh2. 

Shi = ^2Rc{{KiLL + KlRB.)-iKiLR + KiRL)}, 
Sh2 = ^2Re{{K2LL + K2RR)-{K2LR + K2RL)}- 

(63) 

Therefore an evaluation of Ki2pp' determines Sh. 
First we discuss Kippi 

Kipp'=-t J2 ^--'Vp'^G^'Ukm%-;,,io), 

kk' aa' 
^-"-Yl ^•^•^j -(^^Pk^p't'pki^l)Gpk'i^2)^ 

(64) 



kk' 



where duji2 — duJidijJ2. 

Now we examine the fc-sum in ea. (|64p . Appling eq. 
for the summation over /c', essential part of calculations 
are reduced to a form 

EvG;'tpfeK)GjpT,,(^2) 

k' 

- E iVfcf .9X'.'(-i)5p4'fc'(-2)- (65) 

k' 

Notice that the analytic properties of gpf/^p^ki are deter- 
mined by ±ie because of the free electrons in leads. Sum- 
mation of the product of two g° with different frequencies 
LUi and LU2 over k' gives rise to a J-function singularity 
when LUi is close to UJ2, and a principal integration. This 
point is quite different from the one in Si . 

In contrast, the summation over k in eq. (|64p has the 
same character as the one previously discussed in Si. 
Thus it does not produce the J-function singularity. 

To determine Shi we evaluate 2 Re{KiLL + Kirr) 

2 Re{KiLL + Kirr) = 

2 / ^ {T(u;)(l - Tiu;)){fL - JrY + a{u:) - fo(c.)} . 



27r 



(66) 



The term of r(l — T) and the a(a;)-term originate form 
the (5-function term where 

a{u:) = ^2p(^)2(f i - f «)-ii^(/i - !r). (67) 



The &(a;)-term originates from the principal integration. 
We can prove that &(w) completely agrees with a{tj}): 
b{u}) — a{uj). As a consequence 2 Re(i^iii -I- Kirr) is 
characterized by ^^(l — T), as expected for shot noise. 

Moreover the same type calculations lead to that 
2 Re{KiLR + Kirl) - -2 Re{KiLL + Kirr). Shi be- 
comes 



4e^ 

Jhl — —r- 
h 



dujf{L0){l~f{uj)){fL-.fRf 



(68) 



We conclude that Shi has a shot-noise form, but char- 
acterized by the full T. Employing eq.([33]), the leading- 
order of Shi can be calculated at zero temperature as 



where U/'n:T = R-l. 

Having addressed Shi for the bubble-diagrams, let us 
turn to Sh2 for the vertex-diagrams. We can confirm 
that for Sh2 only the term (a) in FigU] remains in the 
same way as S'2. Considering this point, the second-order 
correction to if 2pp' are evaluated. The obtained K2pp' for 
p,p' — L,R decide Sh2 as follows, 



Sh2 = Y I '^■^i'^'^2 ReFh. 



(70) 
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Fh is given by 




p/ 


'^gdp'k'iri-O )gp,k'd{0 -T2) 
k> 


X5d(0+ - n) ^ <P y^ Vpk9dpkiT2 - 0+ 

n h 




X C/gd(ri - T2)Ugd{T2 - n) 



3.4. Fano factor Fh 



In this section, we would like to argue the Fano factor 



Fh 



Sh 
2^ 



(75) 



(71) 

FoUowing the discussion of the fc-sum in Shi, we find 
that [J2k' •■•] ™ eq. ([7T|l gives the (5- function part and the 
principal- integration part. 

Taking this point into account, we perform calcula- 
tions of all Keldysh components of Ff , defined by 



Conventionally, noise power S at zero temperature is 
treated as shot noise. However, we have proposed Sh 
as shot noise at any temperature based on the nonequi- 
librium Kubo formula. Therefore, it is natural to define 
the Fano factor with St in eq. ([75]) . Here h expresses the 
backscattering curren t^'^i^'^'^^i^^ 



2e^ 



(76) 



Ti = t" and T2 = ^2 ™ 6q. ([TT]) • We discuss rele- The / in eq. (gS]) gives an expression of 2e/h 



for Sh2- We find that the principal- 

7^0 



vant Ea/jRei^."^ 

integration parts in X]a/3 ^^^h'^ vanish. Concerning the 
(5-function part, ReF^^^ , ^gF^ and ReF^"*" do not con- 
tribute to the leading-order of Sh2 in V . ReF^"*" deter- 
mines the leading-order of Sh2 

4e2 
^9d 



S, 



h2 



dujiduj2{fL{'^i + UJ2) - fniuJi + cc)2)) 
wi)(/l(wi) - /fl(wi))g''(wi) X r(w2). 



(72) 

where c' = 7rr^C7^/(27r)^. r(a;) is the vertex function 
defined in eq. (|58p . Employing the asymptotic form of 
r(w), Sh2 is evaluated as 

2^ 



4e" 



1 



'/i2 



\v\<^{R-iy 




(73) 



where U/ttT = R-1. 

Therefore, the asymptotic form of Sh = Shi + Sh2 is 
determined as 



Sh 




1) 



eV_ 

T 



(74) 



Finally we summarize the results of the noise power 
S — S1+S2 and the shot noise Sh ~ Shi+Sh2- Subscripts 
of 1 and 2 represent the contributions from the bubble- 
diagrams and the vertex-diagrams. 

As a check we begin with S — Si + S2- Taking the s-d 
limit leads to F ^ Tr- and the Wilson ratio R ^ 2. We 
find that the resulting asymptotic form of 5 = 6*1 -I- S'2 
precisely agrees with the ones obtained as shot noise by 
using the fixed-point Hamiltonian^^i^iiS^. Up to 0{V^), 
the second-order RPT indeed gives the correct result. 

We turn to Sh proposed as the new formula of shot 
noise. Up to the leading order in V, Si = Sih and 6*2 = 
S2h are clearly satisfied. Therefore for the symmetric 
case at zero temperature, Sh precisely agrees with S up 
to 0{V^). 



4e3 I 1 feVV 5 , ,, feVV] , , 

2^^^ = -^|t2(t) + T2(^-i) (t)|-(^^) 

The backscattering current 2el\, just corresponds to Sh\ 
for the bubble-diagrams in eq. (f69|) : 2e/fc = Shi for F > 0. 
Consequently Ff, = {Shi + Sh2)l^e.h = l + S'h2/2e/f,. Us- 
ing expressions of £'^2 for the vertex diagrams in eq. ([75]) 
and 2elh in eq. ((77)) . we obtain F\, 



Fh = l 



^{R-\f 
H-5(i?-l)2' 



(78) 



The first term in Ffc originates from 5^1. Shi is char- 
acterized by only T . If correlation effect gives only a 
change from the bare T" to the full T, F(, = 1 would hold 
within Oiy^^. Actually, Sh2 from the vertex-diagrams 
contributes to F(,. As a result it gives an enhancement 
factor of the second term in eg. ([75)1 . The vertex-diagram 
is known to describe a kind of the back-flow effect. The 
back-flow effect enhances the Fano factor for backscat- 
tering current. The form of the enhancement in Ff, has a 
universal feature that only the Wilson ratio determines 
it for any U . 

Here we check the resulting F(, in the limit of [/ ^ 00 
and C/ ^ 0. Appling limiting values of the Wilson ratio 
R for Fh leads to 



Fi 



5 
) 3 



C/^00 i? = 2 
L/^0 R=\ 



(79) 



Ff, indeed reproduces the universal fractional value of 
Ff, = 5/3 derived in the s-d limit^^ J^^'^^ , and a naively ex- 
pected value Ffc = 1 for a noninteracting system. There- 
fore Ff, given by eq. (j78p is an extension for any I] from 
the already obtained value of F;, — 5/3 in the s-d limit. 

Before closing discussion, we touch on an application 
of Ff,. By using Ff,, we may determine the Wilson ratio 
directly from experiments as follows. 



R= 1 



fb-1 
9 - 5Ff, ■ 
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4. SUMMARY 

The noncquilibrium Kubo formula Sh = S — Ak^TG 
allows us to propose Sh as the new definition of shot 
noise in general. Experimentally measured S — Ak^TG 
can be compared with a theoretical prediction of Sh at 
any temperature. Therefore, the noncquilibrium Kubo 
formula thus opens a new approach to studies of shot 
noise in correlated systems at any temperature and any 
bias voltage. 

Then, using this approach shot noise through a quan- 
tum dot in the Kondo regime has been investigated. For 
simplicity, the symmetric case has been discussed at zero 
temperature. At T = 0, the noncquilibrium Kubo for- 
mula gives Sh = S. We have thus analyzed both S 
and Sh, up to 0{V^) which is the leading order in V 
for the symmetric case. The renormalized perturbation 
theory (RPT) has enabled us to obtain the exact asymp- 
totic form of S and Sh- Both of S and Sh are expressed 
by two-particle Green functions. Thus they have two 
types of contributions from the bubble-diagrams and the 
vertex-diagrams as S* = 5*1 -I- 52 and Sh — Shi + Sh2- It 
has been shown 5i — Shi and ^2 = 5^2 up to 0{V^). 
We have concluded that Sh indeed equals S* at T = 



which was conventionally defined as shot noise. Finally 
we have pointed out that Sh2{— S2) from the vertex- 
diagrams leads to the universal enhancement factor in 
the Fano factor: Ft = l + A{R-l)'^/{l + 5{R-lf). This 
expression includes the result of Ff, = 5/3 in the Kondo 
limit by using the Wilson ratio R — 2. 

Furthermore, we have found that 5*1 splits into the 
thermal-noise part and the shot-noise part with the full 
transmission probability T. Here, let us recall the fact of 
S = StiT) + SsiT) + AS. In fact. Si = St{T) + Ss{T). 
Thus, AS describes nothing but the contribution from 
the vertex-diagrams: ^2 = AS'. As discussed previously, 
concerning AS it is not clear whether noise power S" at a 
general temperature is split into the thermal-noise part 
and the shot-noise part. However, it is impossible to 
ignore the effect of AS, even at zero temperature because 
Sh2 — S2 = AS is always relevant as in the universal 
enhancement factor in Fh. 
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